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ABSTRACT 

We analyze the asymmetry in the partial widths for the decays 5 ± — ► 
MM7r ± (M = 7r + , K + , 7T°, r)), which results from the interference of the non- 
resonant decay amplitude with the resonant amplitude for — > Xco^ fol- 
lowed by the decay Xco — ► MM. The CP violating phase 7 can be extracted 
from the measured asymmetry. We find that the partial width asymmetry for 
B^ — > 7i + 7i^7i ± is about 0.33 sinj, and about 0.45 siwy for B ± — > K + K^ti ± , 
while it is somewhat smaller for B ± — > 7r 7r°7r ± and B — > 777771" . Potential 
sources of uncertainties in these results, primarily coming from poorly known 
input parameters, are discussed. 



The measurement of CP asymmetries in charged B meson decays might 
provide the first demonstration of CP violation outside the K system [I], §]. 
Among the usual three CP-odd phases a, (3 and 7, the phase 7 seems to be the 
most difficult to explore experimentally |J. One possibility to measure this 
CP odd phase 7 = arg(V* b ) has been suggested in j|, |j. In this approach 
the asymmetry appears as a result of the interference of the nonresonant 
— ► MM7r ± decay amplitude and the resonant -B 1 * 1 — > Xco^ ~ > MMtt 1 * 1 
amplitude, where Xco is the ++ cc state at 3.415 GeV. The absorptive phase 
necessary to observe CP violation in the partial width asymmetry is provided 
by the Xco width. There have also been suggestions to determine the CP 
violating phase 7 by analyzing the Dalitz plots in the decays B — > 7i + 7i^7i ± 
§ and B ± -> 7r-7r+is: ± 0. 

On the experimental side the CLEO collaboration has reported upper 
limits on some of the nonresonant decays of the type B + — > h + h + h~ [R]. 
CLEO found the upper limits on the branching ratios BR(B + — > 7r + 7r~7r + ) < 
4.1 x 10~ 5 and -> K + K~t: + ) < 7.5 x 1(T 5 . In @] the branching ratio 

for B + — > 7T + 7r _ 7r + was estimated to be in the range 1.5 x 10~ 5 to 8.4 x 10~ 5 . 

Motivated by this theoretical expectation and the CLEO experimental 
results, we further investigate the asymmetry in the decays 5 ± — > MMir , 
where M = tt + , K + , tt°, 77, improving upon the calculation of the nonresonant 
part of the decay amplitude. We will assume, as in the resonant decay 
amplitude in B ± — > MMtt^ is due to the cc resonance Xco which subsequently 
decays into MM, where M = tt + , K + , 7r°, 77. Note that we are interested only 
in the kinematical region where the MM invariant mass is close to the Xco 
mass, as in [|J. Thus MM arising from other resonances such as the p need 
not be considered. However, we will use the nonresonant B^ — > MMtt^ 
decay amplitudes, calculated using techniques developed previously in our 
analysis of Dn decays |I0| . In particular, we use the factorization approxi- 
mation, in which the main contribution to the nonresonant 5 ± — > MMk^ 
amplitude comes from the product < M M\{ub)y -a\B~ > < n~\(du)v-A\Q > 
where [q^q-^v-A denotes (?i7 /t (l — 75)92- F° r the calculation of the matrix 
element < MM\(ub)v~A\B~ > we extend the results obtained in [|10| , where 
the nonresonant D + — > K~7t + lv decay was analyzed. In this analysis the ex- 
perimental result for the branching ratio of the nonresonant D + — > K~tt + Iv 
decay was successfully reproduced within a hybrid framework which com- 
bines the heavy quark effective theory (HQET) and the chiral Lagrangian 
(CHPT) approach. 
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The combination of heavy quark symmetry and chiral symmetry has been 



quite successful in the analysis of D meson semileptonic decays [ I I ] - [16 



The heavy quark symmetry is expected to be even better for the heavier B 
mesons |L4], [L5j . However, CHPT could be worse in B decays due to the large 
energies of light mesons in the final state. It is really only known that the 
combination of HQET and CHPT is valid at small recoil momentum. In [[16 



we have modified the hybrid model of [3 f| - |14| to describe the semileptonic 



decays of D mesons to one light vector or pseudoscalar meson sate. Our 
modification is quite straightforward: we retain the usual HQET Feynman 
rules for the vertices near and outside the zero- recoil region, as in |14|, 115 ], 
but we include the complete propagators instead of using the usual HQET 
propagator [|16|| . This reasonable modification of the hybrid model enabled 
us to use it successfully over the entire kinematic region of the D meson weak 
decays |TI], [T^]. The details of this approach can be found in [|TI], [Tfl] . 

In the following we systematically use this model [K], |16| to calculate the 
nonresonant 5 ± — > MMw^ decay amplitude. We find there are important 
contributions, which were not taken into account previously H. 



The weak effective Lagrangian for the nonleptonic Cabibbo suppressed B 
meson decays is given by [§j 

£» = -^OWa^'O! + a e 2 ff 2 ) (1) 

where 0\ = (ub)v-A(du)v^A an d 2 = (uu)v-A(db)v-A- We take the ef- 
fective coefficients a\ * (i = 1,2) from the phenomenological fit [[TJ|, which 
gives a\ ^ ~ 1.08 and a e 2 ^ — 0.21. These values are also in agreement with 
other analyses [|l^, p0[| . We do not take into account the contributions aris- 
ing from penguin operators, since these contributions are not expected to be 
important j|, |], [Hj]. The quark currents required in the weak Lagrangian 
([D can be expressed in terms of the meson fields, as previously described 
explicitly in |]| 0. 

Using the factorization approximation, as in Q, we can analyze all possi- 
ble contributions to the MMn^ nonresonant amplitude. The various 
kinds of contributions are shown in Fig. 1. To illustrate the use of the fac- 
torization of the amplitude we consider the specific decay B~ — > 7r~7r~7r + . It 
is easy to see, then, that the two contributions shown in Fig. la, which come 
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from the operator 0\, cancel each other in the chiral limit — *> 0, since in 
this limit 

< ■K~(p 1 )7r + (p 2 )7r'(p 3 )\Oi\B-(p B ) >=< ^(pi)7v + (p 2 )n-(p 3 )\(du) A \0 > x 

<0\(ub) A \B-(p B ) > + 

< TT-(p 1 )TT + (p 2 )TT-(p 3 )\C s \-K-(p B ) > — X 

m B — 

<n-(p B )\(du) A \0><0\(ub) A \B-(p B )> = 

(2) 

^From Fig. lb there is a contribution from operator 2 : 

<7i~(pi)7r + (p 2 )Tr~(p3)\0 2 \B-(p B ) > = 
< 7i~(pi)7i + {p 2 )\{uu) v \0 >< n-(p 3 )\{db) v \B-(p B ) > + 

(Pi^Ps)- (3) 

The matrix element < n~ (pi)n + (p 2 )\(uu)v\0 > can be calculated using the 
model developed previously |TjJ. In this model the matrix element ([J) is 
determined by the pion form factor, which is dominated by the p meson pole. 
However, the p meson contribution is not relevant in the present calculation of 
the nonresonant decay amplitude. Therefore, the contributions coming from 
diagrams in Fig. lb can be neglected in our calculation of the nonresonant 
decay amplitude. 

The only important contributions to the nonresonant decay amplitude 
comes from the diagrams in Fig. lc and is given by 

< n-(p 1 )n + (p 2 )n-(p 3 )\0 1 \B-(p B ) > = 
< n-(p 3 )\(du) V - A \0 >< iC Xpi)it + (p 2 )\(ub) V - A \B- (p B ) > + 

(Pi^Ps)- (4) 

For the matrix element of < K~(pi)it + (p 2 )\(ub) v ^ A \B~{j) B ) > we use the 
results obtained in the analysis of the nonresonant D + — > tt + K~Ii>i de- 
cay width Jl0[. Following this analysis we write the matrix element < 
7r~ (pi)n + (p 2 )\{ub)v- A \B~ (p B ) > in the general form 

< n-(p 1 )n + (p 2 )\u'j IM (l - ^)b\B~(p B ) > = ir(p B - p 2 

+iw+(p 2 + + iw-{p 2 - - ^he m ^p a B p 2 p\ . (5) 
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In the present case only the nonresonant form factors w™' and contribute. 
Note that the contribution proportional to r is of order m\ and therefore can 
be safely neglected. However, in |J, where a different parametrization of 
the form factors was used, in neglecting the contributions of the order 
the contributions proportional to w± were also dropped. In the notation 
of § the product < M(p 1 )M(p 2 )\{ub) v ^ A \B(p B ) > < n(p 3 )\(du) v „ A \0 > is 
proportional to to F^m^ |22j and one can easily show that F^m^ = —m%r + 
ml/ml(w + +w^)(p B +P2)-P3+ m l/ m l( w +— W-)(Pb+Pi)-P3, which explicitly 
demonstrates that the terms proportional to the form factors w± arising from 
the product < MM\(ub) v _ A \B > < iv\(du)v-A\Q > can not be neglected, but 
are important contributions to the nonresonant decay amplitude. Moreover, 
this contribution cannot even be treated as being constant |6], ^] ; it depends 
significantly on the variables s = (pb — P3) 2 = (P2 + Pi) 2 , t = (pb — Pi) 2 = 
(P2 + P3) 2 and u = {p B - p 2 ) 2 = (pi + P3) 2 . 

Using the preceding analysis we can write the amplitude for the nonres- 
onant decay B~ — ► 7r + 7r~7r _ 

M nr (B~(p B )^ir-( Pl )n-(p 3 )ir + (p 2 )) = ^=V* d V ub 

{at ff [^(ml-s-ml)wl r ( Sl t) + ^(2t + s-ml - 3m>!T(f)] 

+ (*«->*)}, (6) 



where 



, 3/2 1/2 1 

nri ,\ 9 jB*7TT'B* m B n I 2 2 j.M 

w + (s,t) = - — [i--^(m B - m ^- t )] 

Jw 1 ~ m B* zm B* 

+ ^-^^-(a + .-rfB-*^), (7) 

Z Jir Z Jtt m B 



and 



n 3/2 1/2 1 , 

T 2 t-m%, [1 + 2<; (m * ^ t)J + _ ^^ 



The parameters 0^2 are explicitly defined in [[T(|. Note that both the «i and 
«2 terms are important in (|7|) and (§), which was overlooked previously ||. 

For the nonresonant decay amplitude B~(ps) — > 7r°(pi)7r°(p2)7r _ (p3) there 
are two contributions: one proportional to a\ < it~\{du)j$ > < 7r 7r°|(u&)^|£>~ > 
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and another one proportional to afe < k°\(uu)a\0 > < 7r°7r \{db)A\B >. 
The amplitude for B~(p B ) — > 7i (pi)7i (p 2 )iT~ (P3) is then given by 

M nr (B-(p B )^ir ( Pl )ir (p 2 )7r-(p 3 )) = ^V ub V^ 

{a[ ff [w n + r (s,t){m 2 B -ml ~ s) + w m \t)(2t + s - m 2 B - 3ml)} 
+a e 2 ff [w r l r (u,t)(m 2 B -ml-u)+w nr (t)(2t + u-m 2 B - 3ml)} 

+ (t^u)}. (9) 

The form factors w± appearing in the part of amplitude proportional to a e 2 ^ 
are given by (^) and (§), including the terms proportional to a 1)2 , while in 
the part of amplitude proportional to a\ these terms depending on at\ j2 are 
absent. 

A similar analysis of the nonresonant amplitude for the decay B~(p B ) — > 
K + (pi)K-(p 2 ) n-(p 3 ) gives 

M nr {B-(p B ) - K+( Pl )K-{p 2 )7r-(p 3 )) = ^V: d V ub 

{at ff [^(m 2 B ~s-ml)wl r (s,t) + ^(2t + s-m% 

-ml - 2m 2 K )w™(t)}}, (10) 



where 

f 3/2 1/2 , 

nrf ,\ 9 JB* m B* m B r-| 1 / 2 2 j.M 

w + (s,t) = - — [l--^(m B -m n -t)} 

Ik 1 ~ m B* zm B* 



+ h-^T 1 ^^- »'L-»r;-2,„7j. (ID 

and 



w nr(t) ~ 9 f B * m B* m B 2 {1 , 1 (n) 2 2 f)] V^B®1 f i n"\ 

The analysis of the decay B~{p B ) — > ^(pi)^^) 71 " ^) is a little more 
complicated due to r] — rf mixing. The nonresonant decay amplitude is 



V2 8 /»? Jr, 



2 
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{a e 1 ff [wl r (s,t)(m 2 B -ml-s)+w n _ r {t){2t + s-m 2 B -ml - 2m 2 ,)] 
+a e 2 ff [wl r (u, t)(m 2 B -m 2 v -u) + w m (t)(2t + w - m| - 2m 2 - m 2 )] 

+ (t <-> u)}.(13) 

The form factors lo^f i n (0) are given by equations (|i~l~l) and (1^), with 
m# — > and without the terms proportional to ax,2- The 77 — rf mixing is 
defined as usual with rj = rj 8 cg — r] sg and rf = r] 8 sg + r] cg, where eg = cosO 
and sg = sinO and 9 ~ —20° [Ell. In the numerical calculations below we 



used the values f v = 0.13 GeV and f v > = 0.11 GeV determined in |p4j . We 
will not analyze cases where the the 7/ meson is in the final state since one 
expects that in the nonleptonic decays of B mesons into final states rfX the 
gluonic penguin contributions are probably very important (see, for example 

0- El)- 



The partial width for the nonresonant decay B — ► MMix is 

T nr (B- -> MMvr-) = / IMuf ds dt. (14) 

(There is an additional factor of 1/2 for the decays with two identical pseu- 
doscalar mesons in the final state.) The lower and the upper bounds on s 
are s min = (m 2 + m 3 ) 2 , s max = [ra B — mi) 2 , while for t they are given by 

t m in,max{s) = ml + m 2 - -[(s - + to 2 ) (s + to 2 - ml) 

T Ai(s,m|,m?)A3(s,77^,ml)], (15) 
where A(a, b, c) = a 2 + b 2 + c 2 — 2(ab + ac + ab). 



Unfortunately there is no experimental measurement of the B meson de- 
cay constant and therefore we will use the usual heavy quark symmetry 
relation [jFJI 

JD V m B 

to obtain fs from the fo, even though it has not been experimentally mea- 
sured either. But, there are some data on fo s @, albeit with large uncer- 
tainty, and taking fp ~ 200 MeV is reasonable |I7|, |2"7fl . Then the B decay 
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the D* 
g = 0.15 ±0.08 

' 9 



constant is fs — 128 MeV. We will use the value of \V u b\ = 0.0031, the 
latest average in ||. The parameter g in the form factors, determined from 
Dtt decay, is g — 0.3 ±0.1 [28|. /,From D° — > K~l + v we found 
[17| . In the present calculations we will consider the range 
0-2 < g < 0.23, the overlap between these two determinations of g. We 
have previously determined the numerical values of cii and «2 analyzing the 
D 



Vlvi decays [16 



|, and must extrapolate to B mesons from these D 
meson values. To this end we apply the soft scaling of the axial form factors 
A\ and A2 as discussed in |14| . This scaling procedure has the virtue that it 
does not neglect the masses of the light vector mesons when the heavy quark 
symmetry is used. It is completely in the same spirit as the basic assump- 
tion underlying our simple modification of the HQET propagators developed 
previously |L0|, [16], [L7j and used in the present analysis. Soft scaling |14| of 



the axial form factors means that 



2 ) 

max J 



const 



II 



and 



A 2 V (Qmax) = COnst 



TUtf + my 



H 



(17) 



(18) 



where H and V denote heavy pseudoscalar and light vector mesons, respec- 
tively. It is easy to see that this scaling leads to the following relations: 
af p and a^* /m D = a^ p jm B . Among all cases found in [16|, [17 



a 



1 



since 

+^ < 



we select the values af K * = -0.13 GeV 1 / 2 , a^ K * = -0.13 GeV 1 / 2 
only this choice of parameters gives 3.4 x 10~ 5 < BR(B~ — > 7r~7r + 7r _ 
3.8 x 10~ 5 , consistent with recent data ||. All the other combinations of 
ai,2, found in |l(J give a B~ — > 7r _ 7r + 7r + branching ratio larger than the 

|. And for this same set of parameters, we also 
^ K~K + n~) < 1.5 x 10~ 5 , which is bellow the 
J. We also note the following limits for the un- 
1.5 x lO" 5 < BR(B- -> n-TT°ir°) < 1.7 x 10~ 5 
-> 7i^r]i]) < 1.1 x 10~" 5 . We note that the contri- 



] give a B~ 
experimental upper limit ||. 
find 1.4 x 10- 5 <BR(B- -> 
experimental upper limit M. 
measured branching ratios: 



and 1.0 x lO" 5 < BR(B- 
butions to the branching ratios arising from a± t 2 are very important in these 
numerical results. 

In addition to the uncertainties in our results arising from the uncertain- 
ties in the values of the parameters discussed above, there is a potentially 
quite large error that could come from the uncertainty in the CKM matrix 
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element V u b and the decay constant fs- For example, the range of values 
0.0018 < V u b < 0.0044 || could change the branching ratios by as much as 
a factor of 2, but the resulting uncertainty in the CP asymmetry, which we 
discuss next, is somewhat smaller. 



In order to obtain the CP violating asymmetry, one also needs to calculate 
the resonant decay amplitude B~ — > XocK~ — * MMtv~ . This amplitude can 
easily be determined in the narrow width approximation, as in ||]: 

M r (B~ -> Xocir~ -> MMtv") = 
M(B- - xocvr-) 1—^= M( Xoc - MM) + (s <- t). (19) 

In our numerical calculations we will use the estimate BR(B ± — > Xco 7r± )/ 
BR(Xco 7r + 7r~) = 5 x 10 -7 derived in ||. The Xco decay data |{J then fix 
the decay amplitudes for Xco MM, (M = n + , ir°, K + , rf). 



Finally we can calculate the partial width asymmetry in the B^ — > 
MM7r ± decays. We are only interested in the kinematical region where 
the MM invariant mass is close to the mass of the Xco meson, m XOc = 



3.415 GeV. The partial decay width T p for B~ 



MMtx , which con- 



tains both the nonresonant and resonant contributions, is obtained then by 



integration from s r 



[m 



A (ic. 



- 2T 



XOc 



to s. 



(m 



AO. 



+ 2T 



A'Oc 



XOc 



0.014 ± 0.005 GeV is the width of the xoc- 



(27r) 3 32m| 



ds 



dt \M nr + M t 



where 



(20) 



Similarly, Tp, the partial decay width for B + 
the nonresonant and resonant contributions, 
is defined by 

A 



r — r- 

I V Pi 

' r +r- 



— > MMtt + , also contains both 
The CP-violating asymmetry 



(21) 



For the range of values of g and selected a.\^ discussed above we obtain 
the ranges 



0.33 sin-f < A(B Z 



TV TT TV 



< 0.34 siwy, 



(22) 



8 



0.44 sirvy < A{B ± -> K + K^tt ± ) < 0.45 sin-/, (23) 

0.23 sirvy < A{B ± -> 7r°7r°7r ± ) < 0.24 sin-/, (24) 

and 

0.17 sinj < A(B ± -> 777771^) < 0.20 sm7. (25) 

In [g it was found that A(B ± -> 7r + 7r"7r ± ) = (0.44-0.49)sm7, which differs 
from ( P2|) due to the importance of the a\£ terms. 

The uncertainties due to the experimental errors in the remaining input 
parameters have not been included here, but we can roughly estimate that 
the rather large current uncertainties in V u b, T Xc0 and T(B~ — > Xco^) could 
result in the error in the asymmetry being as large as even 100%. 

To summarize, we have analyzed the partial width asymmetry in — > 
MM7r ± decays (M = tt + , K + , 7r°, rj), which signals CP violation, and can 
potentially be used to determine sinj. The asymmetry results from the 
interference of the nonresonant decay amplitude with the resonant decay 
amplitude B^ — > Xoc^ followed by xoc MM. The asymmetry, which is 
rather sensitive to the choice of parameters, was estimated to be 0.33 sin 7 
for 5 ± — > 7T + 7r _ 7r ± and 0.45 sin 7 for B^ — > K'K + tt ± , while it is smaller 
for B ± -> 7T 7T 7r ± and B — > t]7]tc decays. The estimates of these partial 
width asymmetries, while perhaps uncertain by as much as a factor of 2, do 
provide useful guidance for the experimental searches for CP violation and a 
measurement of the phase 7. 

This work was supported in part by the Ministry of Science and Tech- 
nology of the Republic of Slovenia (B.B., S.F., and S.P.), and by the U.S. 
Department of Energy, Division of High Energy Physics under grant No. 
DE-FG02-91-ER4086 (R.J.O.). S.F. thanks the Department of Physics and 
Astronomy at Northwestern University for warm hospitality. 
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Figure Caption 

Fig. 1. Skeleton diagrams for the various contributions to the nonres- 
onant B~ — > MMn~ amplitude. The square in each diagram denotes the 
weak transition due to the weak Lagrangean C w (|l|), while each dot denotes 
one of the two corresponding weak currents. 
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